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Abstract 

Responses of small open oscillator systems to applied external forces have been studied with the 
use of an exactly solvable classical Caldeira-Leggett (CL) model in which a harmonic oscillator 
(system) is coupled to finite A^-body oscillators (bath) with an identical frequency (a;„ = oJq for 
n = 1 to N). We have derived exact expressions for positions, momenta and energy of the system 
in nonequilibrium states and for work performed by applied forces. Detailed study has been made 
on an analytical method for canonical averages of physical quantities over the initial equilibrium 
state, which is much superior than numerical averages commonly adopted in simulations of small 
systems. The calculated energy of the system which is strongly coupled to finite bath is fluctuating 
but non-dissipative. It has been shown that the Jarzynski equality (JE) is valid in non-dissipative, 
non-ergodic open oscillator systems regardless of the rate of applied ramp force. 
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I. INTRODUCTION 



In the last decade, a significant progress has been made in our understanding of nonequi- 
librium statistics. Experimental and theoretical studies have been developed on small sys- 
tems such as quantum dots and biological molecular machines which generally operate away 
from equilibrium (for reviews, see Refs. [l 3|). The development of modern techniques 
of microscopic manipulation has promoted experimental studies of small systems. It has 
become possible to study the response of small systems to applied external forces. In par- 
allel theorists have developed the important three theorems: the Jarzynski equality^ 
^, the steady-state and transient fluctuation theorems @-7|, and Crook's theorem 
These fluctuation theorems may be applicable to nonequilibrium systems driven far from 
the equilibrium states. In this paper we pay our attention to a remarkable JE given by 
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(1) 



where W denotes a work made in a system when its parameter is changed, the bracket (■) 
expresses the average over the work distribution function (WDF) P{W) of work performed 
by a prescribed protocol, AF stands for the free energy difference between the initial and 
final equilibrium states, and (3 (= l/ksT) is the inverse temperature of the initial state. 
Equation ([T]) includes the second law of thermodynamics: (yV) > AF, where the equality 
holds only for the reversible process. The JE was originally proposed for classical isolated 
system and open system weakly coupled to baths which are described by the Hamiltonian 
jj] and the stochastic models j8|. Jarzynski later proved that the JE is valid for strongly 
coupled open systems {q]. A generalization of the JE to quantum systems has been made in 
Refs. H-Q. 



A validity of the JE has been confirmed by some experiments 18l-l23[|. Liphardt et. 
al. [isl have determined the free energy required to unfold a single RNA chain from non- 
equilibrium work measurements. Wang et. al. [19] have considered a colloidal particle pulled 



20 



2l| have verified the JE for 



through liquid water by an optical trap. Douarche et. al. 
a mechanical oscillator that is driven out of equilibrium by an external force. By using a 
torsion pendulum composed of a brass wire, Joubaud et. al., 



22 



23| have experimentally 



studied the JE of the harmonic oscillator in contact with a thermostat and driven out of 
equilibrium by an external force. 
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Some criticisms, however, have been raised for the vahdity of the JE 2J]- 35| . Cohen and 



Mauzerall j2J] pointed out that it is difficult to define the distribution and the temperature 
duxing the i— le p.oce.. In „ to tMs cHt.cn, Ja..ynsM ha. cWn,ed that 
the JE holds if the initial state is in the equilibrium state with the definite temperature M. 



It has been pointed out that t 
and in a rigid rotator model 



ne 



26 



E may be violated in ideal gas model j25|. 



31 



32 



28 



35| 



29(1 . Therefore it is currently an important issue to 



examine the validity condition of the JE. 



Vlaiw 



|3|,Q, 



studies have been reported for harmonic oscillator systems by both experimental 



20 



23| and theoretical methods 36|-46|. Theoretical ana^ 



oscillators with t 
Langevin model 



le use of the Markovian Langevin model 



20 



23 



yses have been made for 
the non-Markovian 



37H40|. Fokker-Planck equation |4l|, and Hamiltonian model 42h46| . All 



of these studies have shown that the JE holds in isolated and open oscillators, assuming 
dissipative memory kernels or the over-damped models. This assumption seems reasonable 
in the situation under which the relevant experiments 20l-l23| have been performed. Recent 
theoretical studies, however, have demonstrated that the energy dissipation is not realized in 
a small system coupled to finite thermal baths 471, |48| . This is quite different from the case 
of infinite baths in which dissipation is realized. Indeed, it is commonly believed that the 
dissipation is realized only when the system is coupled to infinite bath (except for chaotic 
baths) 491]. Poicare recurrence time is finite for finite bath. 

It is necessary to make detailed calculations of responses of small systems to the applied 
force such as variations of position and energy of the system, which have not been reported 
as far as we are ware of. The purpose of the present study is twofold: to make detailed 
study of the response to an applied force and to examine the validity of the JE in open 
harmonic oscillator systems in the non-dissipative situation. We consider the Caldeira- 
Leggett (CL) Hamiltonian model [50|, l5l|, adopting a single-u bath containing uncoupled 
A^-body oscillators with an identical frequency: Un = Uo for n = 1 to N [Eq. ( fT3l) ]. The 
CL model with a single-o; bath is exactly solvable. A similar optic-phonon-mode model for 
bath was adopted in a different context from the present study |52[. In the conventional 
approach, we obtain the Langevin equation from the CL model, with which its properties 
are investigated. In this study, we have directly obtained the Laplace-transformed equation 
of motion of the system. The energy and work of the system induced by the applied force are 
analytically averaged over the canonical distribution of initial equilibrium states. Our non- 
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dissipative system-plus-bath yields non-ergodic solutions, for which the JE will be shown to 



be valid in contrast with Refs. 



39 



40| claiming the importance of the ergodicity. 



The paper is organized as follows. In the next Sec. II, we derive expressions of response of 
positions, momenta and system energy induced by an applied ramp force in open oscillator 
systems, by using the CL model with the single-w bath mentioned above. We obtain the 
WDF and the averaged work with which the validity of the JE have been investigated. Some 
numerical calculations are presented. In Sec. Ill an application of other types of external 
forces to the system is studied. We compare our study with the method using the Langevin 
model derived from the CL model. Sec. IV is devoted to our conclusion. 



II. THE ADOPTED MODEL 



A. Equations of motion 



We consider a system of a classical oscillator coupled to a bath consisting of A^-body 



uncoupled oscillators described by the CL model jSOl. l5l|. 

H = Hs + Hb + Hj, (2) 

with 

- + (4) 

n=l ^ ^ 

(5) 



n=l ^ 



where Hs, Hb and Hj express one-dimensional Hamiltonians of the system, bath and in- 
teraction, respectively, M (m), Q (un), Q (qn) and P {pn) denote mass, frequency, position 
and momentum, respectively, of the system (bath), c„ the interaction between the system 
and bath, and f{t) an applied external force. Equations of motion for Q and g„ are given 
by 

^ ' f{t\ (6) 



MQ = -Mn^Q + c„ (qr. 

n=l ^ 



mqn = -mulqn + CnQ. (7) 



Applying the Laplace transformation to Eqs. ([6]) and ([7]), we obtain 



N 



M[s^Q{s) - Q(0) - sQ{0)] = -MQ^ Q{s) - ^ 



n=l 



N 



n=l 



where 



Q{s) 



dt e-''Q{t), 



Q(s) + 5Zc„g„(s) + /(s),(8) 

(9) 

(10) 



and similar expressions for g„(s) and f{s). Solving Eq. ([9]) in terms of q{s) and substituting 
it into Eq. ([8]), we obtain 



Q{s) = G{s) 



N 



Q{0) + sQ{0) + Y, 



c„[g„(0) + sg„(0)] ^ f{s 



M 



(11) 



where the Green's function G{s) is given by 



N 



G{s) = is' + n' + J2 



els' 



n=l 



(12) 



In order to make analytic calculation feasible, we consider a bath containing A^-body 
uncoupled oscillators with an identical frequency Uo and a uniform coupling Cq, as given by 



N 



for n = 1 to A^. 



(13) 
(14) 



We have chosen c„ such that it yields a non-diver gent result in the limit of A^ — )■ oo in Eq. 
( IT2|) (related discussion being given in Sec. III.B) [53.]. With the use of Eqs. ( [T3|) and (fT4|) . 

Q{s) becomes 



Q(s) = G{s) 



Co 



N 



( PnO 



m 



+ sqno] + 



M 



with 



G{s) 



els' 



(15) 



(16) 



where Pq = MQ{0), Qq = Q{0), Pno = f^Qni^) and g^o = Ini^)- Equation ( TTB]) may be 
rewritten as 



G{s) 



(17) 
(18) 
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with 



° Mmujt 



+ (-1 



'Dr. 



(z = l,2), 



Do = {^'-ujlf + 



hi 



aj — a2 
Then Eq. (fT5|) becomes 

2 



+ 



2 > 0, 



hi 



1=1 



N 



/(£) 
M 



whose inverse Laplace transformation yields 



N N 

Q{t) = ^{t)+XQ{t)Q^ + Xp{t)P^ + Y,{t)Y,(ln^ + yp{t)Y,Pr.O, 

n=l 



n=l 



with 



2 

XQ{t) = hi cos g^t, 

^ _ ^ / 6i Co ^ (cos Wot - COS ait) 
i=i 



MVN 



6j Co \ (oj sin cjo^ — ti^o sin ait) 



mUoai{af - u"^) 



B. Position, momentum and system energy 



(19) 
(20) 
(21) 



(22) 



(23) 



(24) 
(25) 
(26) 
(27) 
(28) 



It is necessary to evaluate physical quantities averaged over the canonical distribution of 
initial states, Qo, Pq, {g„o} and {pno}, of the equilibrium coupled system-and-bath H(t = 0). 
In order to make such evaluations, we need following (fluctuation-dissipation) relations for 
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/(O) = given by 



Mn^Ql), = ^-^ = kBT = ^, (29) 
mijjl{qnoqm)o = keT 5ne + -^^^^-^^r^, (30) 



m 



= ksTSni, (31) 



mu 



{Po Qo)o = {Pq qno)o = {Pq Pno)o = {PnO *o)o = {PnO Q)o = 0, (33) 
(Qo)o = {Po)o = (gno)o = (Pno)o = 0, (34) 



with 



_ Tr {e-mo) 0} 

where O denotes an operator and Tr the trace over initial state of i^(0) with Qq, Pq, {qno} 
and {pno}- Equations (l30ll and ( l32l) arise from the relation, 

- 3% ) ( *o - J \ = ksT 5ni. (36) 
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In the limit of Co = 0, Eqs. f l29|) -f l34|) reduce to the well-known result for isolated system 
and bath. With the use of Eqs. (12^ and (El]), the averaged position and momentum of the 
system are given by 

Q{t) = (Q(t))o = m = f xit - t')f{t') dt', (37) 

P{t) = (P(t))o = Mm = M [ x{t- t')f(t') dt', (38) 

Jo 

with the time-dependent susceptibility x(^)? 

, , bi sin ait . . 

XW = E^J^. (39) 

where dot (■) stands for the derivative with respect to time. It is easy to see from Eq. flTSj) 
that the Laplace-transformed susceptibility is given by x{s) = G{s)/M. The frequency- 
dependent susceptibility xi!^) is given by 



oo 

iuit 



X{uj) = / e-'^'x{t)dt = x{-iuj), (40) 







whose imaginary part becomes 



Im = ^ ( 



TiMbi 
2a,; 



[S{uj - ai) - 5{u + ai)]. 



(42) 



The system energy Es averaged over the initial state is given by 



Mm.\ 



E. 



{Q'),-fit) (g)o. 



(43) 



By using (l29i) -( l34l) and Eqs. ( l43l) . we obtain Es given by 



(44) 



with 



E 



(0) 



+ 
+ 

+ 



2MQ? 
MksT 

2 



MXp{tf + Mn^Xp{tf 



+ 



NmkBT 



MYp{tf + Mfi2rp(t)^ 



1 + 



MXQ{t)Y,{t) + Mn2xQ(t)y;(t) 



(/) 



(45) 
(46) 



Here E^^'^ expresses the system energy depending on the temperature but independent of the 
apphed force: E^'p denotes the response to the force: $(t), -^q(^), Xp{t)^ Yq{t) and Yp{t) 
are given by Eqs. (12^-( 128|) : Xqif), Xp(t), Yq(t) and Yp(t) are their derivatives with respect 
to time. It is noted that Q{t) and P{t) are independent of because of the A^-independent 
$(t) in Eq. (!24|) . Furthermore Es does not depend on because the factor in the fourth 
term of Eq. fHSj) is canceled out by the 1/A^ term in Yq(t)'^ in Eq. fl27|) and because the \/N 
term of the last term of Eq. (I45l) is cancelled out by the l/\fN of Yq{t). These properties 
arise from our adopted model with c„ = Co/\fN in Eq. ( |T4l) jsS]. 

The advantage of expressions given by Eqs. (|57l) . (1551) and is that canonical av- 

erages over the initial state have been analytically made and they are free from the numerical 



averaging which is one of difficulties in direct simulations of small systems [4 



48 



58 



-|60|. 



We have so far not specified the form of an external force f{t). For a while we consider 
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a ramp force given by 



fit) 



for t < 0, 

g (i) for < t < r, 

g for t > T, 



(47) 



where r stands for a duration of the apphed force and g the magnitude of the force at t > r 
For the ramp force, Eq. (1241) leads to 

{ttit — sin ait) for < t < r 



*(*) = E ( 

i=i ^ 



Mafr 



El 



Maf J V a,r 



1 



[ttiT + sin ai(t — r) — sin ajt] J for t > r 



(48) 

(49) 



In the following, we examine the three cases of (1) no couplings (co = 0), (2) transient force 
(r = 0) and (3) quasi-static force (r oo). 

(1) In the case of Co = where Eqs. ffT^ - fl?Il) lead to ai = fl, a2 = ujo, &i = 1 and 62 = 0, 
we obtain 



for < t < r 
[Qt + sin Q{t — r) — sin Qt] 



for t > T. 



(50) 
(51) 



Equations ( 125|) and fl26|) lead to 



XQ(t) = cosi7t. 



sin f2t. 



(52) 
(53) 



Es(t) becomes 

Esit) = ksT 



9' 




2Mfi2 








9' 


1 - 


2Mfi2 





2(1 -cosr^r)" 



for < t < r, 
for t > T. 



(54) 
(55) 



(2) In the case of r = 0, Eq. (09]) yields 

gh 



i=l 



(56) 
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which becomes for c, 



(l-cosQt) for t > 0, 



(57) 



yielding 



Es{t) 



(3) In the case of r — )• oo, Eq. 



kBT 



yields 



for t > 0. 



(58) 



which becomes for Cn 



m - 

0, 



i=l 



gbj 



for < t < oo, 



(59) 
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for < t < cxD, 



(60) 



leading to 



Es{t) 



kBT- 
ksT 



9' 



2 



9 



for t 



r — 7- oo. 



(61) 
(62) 



We have performed numerical calculations for averaged position, momentum and energy 
of the system with M = m = 1.0, Q = Uo = 1.0 and g = 1.0 which are adopted in all 
our calculations otherwise noticed. Position, momentum and energy (work) are measured 
in units of ^JkBT/MVt'^, y/WksT and ksT, respectively. Model calculations of averaged 
positions and momenta are presented in Figs. [l]^a)-(h) where solid and dashed curves express 
Q(t) and Pit), respectively. Figures (U^a) and (b) show the results of Co = 0.0 and Co = 1.0, 
respectively, when a ramp force with r = 100 is applied. Figure [U^a) shows that Q{t) is 
linearly increased at < t < 100.0, and it becomes constant at t > 100.0 where a force g is 
still applied. This behavior is not changed even when the system-bath coupling is introduced 
as shown in Fig. [T](b). Figures [D^c), (e) and (g) show Q{t) and P{t) for ramp forces with 
r = 10.0, 5.0 and 0.0, respectively, applied to uncoupled systems (co = 0.0), where regular 
oscillations are induced. Figures [D^d), (f) and (h), however, show that irregular oscillations 
are induced by external forces with r = 10.0, 5.0 and 0.0 in coupled systems. 

Model calculations of system energy Es(t) are plotted in Figs. [2](a)-(j). Figures [2]^a) and 
(b) show Es{t) for Co = 0.0 and Cq = 1.0, respectively, without external forces [f{t) = g = 0.0] 
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for which Es is constant. Figures |2]^c) and (e) (Figs. EJ^d) and (f)) show Es for Co = 0.0 
(co = 1.0), with apphed forces of r = 100.0 and 10.0, respectively, where Es is gradually 
decreased by an applied force. As far as the uncoupled system is concerned, this behavior 
is not modified when the force with smaller r is applied, as shown by Figs. W^g) and (i) 
for T = 5.0 and r = 0.0, respectively. However, when the ramp force with smaller r is 
applied to coupled systems, the behavior is changed: irregular oscillations are induced in 
Es as shown by Fig. EJ^h) and (j) for r = 5.0 and r = 0.0, respectively. These oscillations 
in coupled systems are realized for ramp forces with t < To where To (= 27i/Q) denotes the 
period of system oscillation. We note in Figs. MM) and (j) that this irregular oscillation is 
not dissipate, which has been confirmed with calculations for t G [0, 10000] (relevant results 
not shown). The averaged system energy in the coupled small systems shows irregular 
non-dissipative oscillations although the total energy of the system-plus-bath is constant 



C. Work and work distribution function 

Next we consider a work performed by an applied external force. By using Q{t) given by 
Eq. (123!) . we obtain the work performed by the force f{t) applied for < t < r ^], 



47|,|48|. 




(63) 



(64) 



n=l n=l 



where 




(for 1] = q and p). 



(for ^ = Q and P), 



(65) 



(66) 



(67) 



With the use of Eqs. (I64D-(I67D, the WDF of P{W) is given by 



P(W) = {5{W-Wo)) 



' 



(68) 
(69) 
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where 



{exp{-iuWo))n = exp(-m0) — 



27r y V 27r 



N 



dQo exp 



- IuCqQq 



X J dPoexp 

N 



- tuCpPo 

2M 



X 



Y[ / dqnoexp 

n=l ^ 

N 



moot, 



iuDgCnQo 



n=l 



(70) 



Performing the Gauss integrals, we obtain 

(exp(-mM/o))o = exp 

where 



—ZU0 — 



u 

2^ 



(71) 



1 



(72) 



With the use of Eqs. (jM]) and ([71]), P(W) is finally given by 



P{W) 



V2^ 



exp 



2^2 



with 



fi = {W) = 0, 



(73) 



(74) 



where and are given by Eqs. f lB^ and fl72|) . respectively. The average of e over 
P(iy) is given by 



which leads to 



fis -iln(e-''"'> 



2 " 2 

It is worthwhile to point out that Eq. ( 1751) may be alternatively obtainable by 



(75) 

(76) 
(77) 

(78) 
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where (f) and are given by Eqs. ( 165|) and ( !72|) . respectively. 
For a ramp force given by Eq. fj47|) . Eqs. fl65|) -f l67|) are given by 



i=l 



1 

2^ 



1^1 — cosoir) 



hi sin ajT 



i=l 

Ml ^ 



COS a,r 



i=l 



Cog 



NM 
Co9 



bi{ai sin — sin ajr) 



NMm 



y-v ^ COSCJoT 

i=l 



cos ajr)] 



(79) 
(80) 
(81) 
(82) 
(83) 



It is noted that R given by Eqs. ([72]), ([77D, (I7S])-([H3]) is independent of because ViV 
factor in the first term of Eq. ( 1721) is cancelled out by 1 / of Dq in Eq. ( I82l) , and because 

factors in the third and fourth terms in Eq. fl72|) are cancelled out by factors of 
D"^ and D"^ in Eqs. (182|) and (!83|) . This is the consequence of our choice of c„ in Eq. (fT5|) : 
a different choice of the dependence of Cn leads to A^-dependent R. Furthermore R is 
independent of (3 because (3 factor in the second term of Eq. ( 1771) is cancelled out by in 
Eq. (Ha). 

Figure[3t^a) shows the r dependence of /i (= (IV)) for Co = 0.0 (sohd curves), 0.5 (dashed 
curves) and 1.0 (chain curves). For r < To (~ 6), we obtain (W) > AF (= —0.5) signaling 
the occurrence of the irreversibility. At the same time, cr (= ^{(W-{W)y)) is rapidly 
increased for r < To, where fluctuation in W much grows, as shown in Fig. [3]^b). For 
r = 2mn/Q {m = 1, 2, ■■) with Co = 0.0, a vanishes [Eq. ( 193|) ]. Figure[3]^c) will be explained 
shortly. 

Figure m shows 3D plots of WDF of P{W) as functions of W and r for Co = 0.0: result 
for Co = 1.0 is not so different from that of Co = 0.0 on first glance. With decreasing r, the 
center of P{W) moves to zero and its width is considerably increased as Figs. [3]^a) and[3t^b) 
show. 
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D. Jarzynski equality 



In this subsection, we consider the JE given by Eq. ([T]). From Eqs. ([T]) and (1751) - (1771) . 
the JE is satisfied if the relation given by 

/? = 0-^ = AF, (84) 

holds. Here AF denotes the free energy difference between the two equilibrium systems with 



and without a force g defined by {q] 



with 



where 



AF = FM-F(0) = -ilu|||. (85) 



U9) = ^ (86) 



n ( ^ P'MQ'f g Y g' 



2^2 



n=l ^ 



TV 



Hb + Hj = Y, 



2m 2 \ " mu^ 



(89) 

Z,lll, z, \ ■iii.uw I 

n=l 

Zs{g) denoting the partition function of the system with Hs{g) for a constant force of 
f{t) = g. By using Eqs. ([85])- ([89]), we obtain 

Zs{9) = (1^) e^^^/--^ (90) 

yielding 

AF = , (91) 

which is independent of the coupling Cq. 

In what follows, we examine /i, cr^ and R in the three limits of (1) no couplings (Cq = 0), 
(2) transient force (r — )► 0) and (3) quasi-static force (r — )► oo). 
(1) In the limit of Cq = 0, we obtain from Eqs. fl72p and fl79p - 



2 2(7^(1 — cosier) 
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leading to 



where AF is given by Eq. (I9T!) . 

(2) In the hmit of r ^ 0, Eqs. ([72]) and ([79])-([83D lead to 

/i = 0, (95) 

where we employ the relations: Cq = —g and Cp = Dq = Dp = 0. A substitution of Eq. 
([96D into Eq. ([77D leads to 

(3) In limit of r — cxd, we obtain 

= -^^'E^ = -^' (98) 
CT^ = 0, (99) 

employing the relations: Yl'i=ii^i / ^1) = 1/^^ ^"^^ = Cp = Dq = Dp = 0. Equations 
dZZD and ([99]) lead to 

Equations ([94]), ([97]), and ffTonl) imply that the JE holds in the three limits of (1) Co = 0, (2) 
r — )■ and (3) r — t- oo. 

Figure [3](c) shows that the JE is numerically verified for 10^^ < t < 10^ with Co = 0.0, 
0.5 and 1.0. The JE is valid even when we adopt other sets of model parameters. It is 
surprising that complicated expressions of fi (= 0) and given by Eqs. ([79]) and ([72]) . 
respectively, satisfy the JE given by Eq. ([8^ . Although the validity of the JE is confirmed 
by numerical calculations, we have not succeeded in its analytical proof except for the three 
cases of Co = 0, r — 7- and r — )► oo. 



15 



III. DISCUSSION 



A. Canonical average over initial equilibrium state 

It should be stressed that the canonical average in Eq. fl5^ must be per formed over the 
total Hamiltonian H (= Hs + Hb + Hj) in the initial equilibrium state [61|. If the average 
in Eq. ( 135|) is performed over the Hamiltonian of the uncoupled state {Hs + Hb) in place 
of if, we obtain a wrong result. Figures [5]^a) and (b) show Esit) (= {Es{t))oo) with no 
forces (/ = 0.0) and a ramp force of r = 100, respectively, with Cq = 1 when the average is 
performed over the initial uncoupled state of {Hs + Hb), 

r^j. {g-/3[Hs(0)+//B(0)] 

^^^00 = Tre-/^I^s(o)+Hs(o)] ' (10^) 

where O stands for an operator. Results in Figs. O^a) and (b) are quite different from 
the corresponding ones averaged over H which have been shown in Figs. E^b) and (d). In 
particular, the irregular energy exchange between the system and bath occurs even when 
f{t) = 0.0 in Fig. M^a), while the initial serene state persists in Figs. EJ^b). Figure [5t^a) 
denotes the result of the case where the system-bath coupling is suddenly added at t = 0.0 
to the uncoupled system in equilibrium state at t < 0.0. The perturbation of the added 
coupling induces the irregular energy exchange between the system and bath which does not 
dissipate for t > 0.0. 

If the canonical average in Eq. fl78|) is performed over Hs + Hb, we obtain 



{e-^^o^oo = e-f'^'^-f'-"l^\ (102) 

with 

^ \ CI „ /VD2 _1 

(103) 



+MCI + ^ + mNDl 



M9? ' moj 



where (f) is given by Eq. f l65|) . Because a' ^ is different from o"^ in Eq. ( !72|) , it wrongly leads 
to a violation of the JE: i? = — fia' 7^ — /3o"^/2 = AF. The related discussion will 
be given also in Sec. III.D. 
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B. A two-step ramp force 



Besides a ramp force given by Eq. ( H7|) . we have employed a two-step ramp force given 



by 



fit) 



\ 'm- 



{l-h)t+{hT-Tm) 
(r — Tm) 



for t < 0, 

for < t < Tm, 
for Tm <t < T, 

for t > T, 



(104) 



where h stands for a magnitude at a middle time of (< r). Figure E] (a)and (b) show fi 
and a, respectively, as a function of r when a two-step ramp input f{t) given by Eq. fll04p 
with g = 1.0, h = 1.5 and = r/2 is applied [/(/!:) is shown in the inset of Fig. EJ^b)]. The 
input force f{t) first linearly increases to 1.5g at t = and then it linearly decreases to the 
final value of (? at t > r. The r dependences of fi and a shown in Fig. are rather different 
from those for a single-step ramp input shown in Fig. El In particular, magnitudes of fi and 
a have resonance-like peaks at r ~ To. Nevertheless the JE holds also for the two-step ramp 
force. 



C. Baths with multiple oj and infinite 

Our study in the preceding section has been made for the CL model with the single-w 
bath, which may be extended to multiple-o; bath. The Green's function given by Eq. ( IT^ 
may be generally expressed by 

1=1 ^ ' ' 

where Ui denotes the normal-mode frequency of the coupled system-plus-bath and 6, is 



62| 



expressed in terms of the corresponding eigenfunction 
The inverse Laplace transformation leads to 

G(t) = X;^^^- (106) 

. -I 
1=1 

The Green's function given by Eq. fllOSp or fll06p has the same structure as that for the 
single-cj bath given by Eq. (|T8|) . Calculations of Q{t), Wq and P{W) may be formally 
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performed in the same way as was made in Sec. II. Then properties of the CL model with 
finite- iV multiple-o; bath are essentially the same as those with single-o; bath. 

On the other hand, in the limit of — )■ oo, the summation over n in the Green's function 
of Eq. f[T^ is converted to integral over a continuous spectrum and it may be expressed by 



G{s) 



s^cl r D(co) 



Mm J w^{s'^ + w^) 



-1 



where D{u) denotes the density of state, 



1 ^ 



(107) 



108) 



n=l 



When we assume the Debye-type density of states: D{uj) = auj^ (a: constant), G{s) is given 
by 

1 



(s + Ci)(s + C2) 



with 



= ±M^^-f^V+f^V 
V V4Mm ) \^Mm J 



(109) 



(110) 



Because of real parts in ci and C2, the inverse Laplace transformation of G{s) in Eq. fll09p 
yields dissipative G{t) which vanishes at t — )■ 00. For dissipation it is necessary that the 
frequencies {un} have a continuous spectrum in the limit of — )■ cxo 62|. With a discrete 
spectrum for finite N, however, the Green's function G{t) in Eq. ( fT2l) is non-dissipative and 
not vanishing in the limit of t — )■ 00. 



D. The generalized Langevin approach 



In the conventional approach to the CL model, we derive the Langevin equation given by 



with 



MQ = -Mn^Q- / ^{t-t')Q{t')dt' + Cit) + fit), 
Jo 



at) 
lit) 



C(t)-7(t)Q(0), 



N 



n=l 
N 



El — ^ ) cosa;„t. 



moot, 



n=l 



g„(0) cos Unt 



PniO) 



mUn 



sin ujnt 



(111) 

(112) 
(113) 

(114) 
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after obtaining a formal solution of g„(0) from Eq. ([7]) and substituting it into Eq. 



50 



5l| . Equations flllip -f lTT4|) express the non-Markovian Langevin equation with colored 



noise. 



When we adopt the single-cj bath given by Eq. f[T5]) . 7(t) and ({t) are given by 

^2 



7(t) 



2 , cos Wot, 



N . , N ' 

EsmWot 
gnO H > ,PnO 



71=1 



mujo 



n=l 



By using the Laplace transformation yielding 

7(s) 



Co 



N 



(115) 
(116) 

(117) 
(118) 



o n=l ■ n=l 

we obtain an equation for Q{s) which is exactly the same as Eqs. f llSp and fll6p . 

It has been shown that the JE is satisfied in the non-Markovian Langevin model with 



colored noise (generalized Langevin model) 



37 



40|, which is different from our Langevin 



equation given by Eqs. fillip . f lll5p and flll6p in the two points: (a) the second term of 
Eq. f lll2p includes an additional term of —'y{t)Q{0) which is missing in the conventional 
generalized Langevin model, and (b) the memory kernel given by Eq. f lll5p is oscillating and 
non-dissipative while that in the generalized Langevin model is dissipative. In the literature 



[e.g. Ref. 



65|), the additional term of —■y{t)Q{0) is discarded and the fluctuation-dissipation 



relation is given by 



(CWC(0)oo = A;BT7(t-t'), 



(119) 



which is derived from the equi-partition relations: (g„o Qeo)oo = {ksT /mu^dni and 
(PnoP^o)oo = {kBT'm)6ne, the bracket (^o denoting the canonical average over the uncoupled 
initial state Hs + Hb [Eq. (llOip ]. If the additional term is included, we obtain 

7(t)7(t') 



(C'W C'(t'))oo 



lit - t') + 



(120) 



which is different from Eq. ( 1119^ . It is noted, however, that when employing Eqs. ( l29l) - (!34ll 
valid for equilibrium initial states of coupled Hamiltonian H (= Hs + Hb + Hj), we obtain 
the desired fluctuation-dissipation relation. 



(C'WC'(O)o = kBT^{t-t'). 



;i2i) 
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Then the ostensible inconsistency of the item (a) may 
importance of the ergodicity is emphasized in Refs. 



3e resolved. As for the item (b), the 



39 



40|] from a study on the JE for 



the generalized Langevin model. Non-ergodic solutions of the non-dissipative generalized 



Langevin equation have been discussed in Refs. 62|-64|. Our non-diffusive memory kernel 
yields non-ergodic solutions for the Langevin equation given by Eqs. fillip . flllSp and f lll6p . 
It is noted that the JE holds in our calculation even if the condition of the ergodicity is not 



satisfied, in contrast with Refs. 



39 



40|. 



Before closing Sec. Ill, it is necessary to mention that Ref. 29|] has studied the validity 



condition of the JE for a general classical dynamical system with any time-dependent exter- 
nal force /(t). It is shown in 29(| that the JE holds for a classical system during a transition 
process in which the value of a parameter / in the system Hamiltonian switches from /q to 
/i in time r, as long as the phase space extension of the system ^q*^ at the initial equilibrium 
phase space with / = /q is the same as the equihbrium phase space Vf^ with f = fi- It is 
noted that the general validity condition reported in Ref. 29|] is satisfied for our system even 



though it is a non-ergodic one. 



IV. CONCLUSION 



We have studied the response to an applied force of small open oscillator system described 
by the exactly solvable CL model with the non-dissipative single-o; bath. Although the model 
adopted in our study seems a pedagogical toy model, it is expected not to be unrealistic 



47 



because non-dissipative properties are realized in small systems coupled to finite baths 
48| . We have obtained exact expressions for position, momentum and energy of the system 
whose canonical averages have been analytically performed over initial equilibrium states. 
Our calculations of system energy and work have shown the following: 

(i) the energy of the system strongly coupled to finite bath is fluctuating but non-dissipative 
in general, and 

(ii) the JE is valid in non-dissipative non-ergodic systems. 



The item (i) supports direct simulations for open systems coupled to finite bath^ 
although it is contrast to the result showing the dissipation for > 10 — 20 



47 



48| 



The 
Our 



item (ii) is consistent with Jarzynski's proof for an arbitrary classical open systems 
study is complementary to the previous studies for dissipative oscillator systems with the 
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use of the Markovian 
equation 4l| and Hamiltonian models 



2Ch23|| and non-Markovian Langevin models 37h40|. Fokker-Planck 



42|-|46|. 



Although the items (i) and (ii) hold for open systems described by the CL 50|, |5l| and 



Ford-Kac models 



65j, it is not certain whether they are valid for any non-dissipative non- 



ergodic open system. In this respect, it would be interesting to examine a work in the 
{Ns + Nb) model for a classical A^^-body system coupled to an A'^-body bath |48|. The 
{Ns + Nb) model clarifies some interesting issues such as the A'^-dependent non-Gaussian 
energy distribution of the system 48] which has been not realized in previous studies for 



CL-type models with Ns = 1 
reported in a separate paper. 



50 



51 



65| . Such a calculation is in progress and will be 
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FIG. 1: (Color online) The time dependence of the averaged position Q (solid curves) and mo- 
mentum P (dashed curves); (a) r = 100.0, Co = 0.0, (b) r = 100.0, Co = 1.0, (c) r = 10.0, Co = 0.0, 
(d) r = 10.0, Co = 1.0, (e) r = 5.0, Co = 0.0, (f) r = 5.0, Co = 1.0, (g) r = 0.0, Co = 0.0, and (h) 
r = 0.0, Co = 1.0. 



FIG. 2: (Color online) The time dependence of the averaged system energy Eg] for no forces 
(/ = 0) with (a) Co = 0.0, (b) Co = 1.0; for the ramp forces with (c) r = 100.0, Co = 0.0, (d) 
r = 100.0, Co = 1.0, (e) r = 10.0, Co = 0.0, (f) r = 10.0, Co = 1.0, (g) r = 5.0, Co = 0.0, (h) r = 5.0, 
Co = 1.0, (i) T = 0.0, Co = 0.0, and (j) r = 0.0, Co = 1.0. 



FIG. 4: (Color online) 3-D plots of P(M/^) as functions of W and r for Co = 0.0, the ordinate of (b) 
being enlarged compared to that of (a). 



FIG. 5: (Color online) The time dependence of Es{t) with (a) no forces (/ = 0) and (b) a ramp 
force of r = 100.0 with Co = 1.0 when the average is performed over initial uncoupled state of 
Hs + Hb- (a) and (b) should be compared to Figs. EJ^b) and (d), respectively, which are averaged 
over the initial coupled state of Hs + Hb + Hj (see text) . 



FIG. 3: (Color online) The r dependence of (a) /i {={W)), (b) a (= {{W - {W))^)) and (c) R 
(= —P~^ln{e~^'^)) for Co = 0.0 (solid curves), 0.5 (dashed curves) and 1.0 (chain curves); arrows 
along the right ordinates in (a) and (c) express AF (= —0.5). In (c) R = AF for Co = 0.0, 0.5 and 
1.0. 
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FIG. 6: (Color online) The r dependence of (a) fi {={W)) and (b) a (= {{W - {W))^)) for a 
two-step ramp force f{t) given by Eq. (|lU4p with g = 1.0, h = 1.5 and = r/2 [see the inset of 
(b)] with Co = 0.0 (solid curves), 0.5 (dashed curves) and 1.0 (chain curves); an arrow along the 
right ordinate in (a) expresses AF (= —0.5). 
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